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Abstract 

An algebra isomorphism between algebras of matrices and difference operators is 
used to investigate the discrete integrable hierarchy. We find local and non-local fami- 
lies of i?-matrix solutions to the modified Yang-Baxter equation. The three i?-theoretic 
Poisson structures and the Suris quadratic bracket are derived. The resulting family of 
bi-Poisson structures include a seminal discrete bi-Poisson structure of Kupershmidt 
at a special value. 
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1 Introduction: A Matrix Formulation of the Lattice Hierarchy 



A popular framework for dealing with infinite lattice systems is the differential-difference 
calculus [I] based on a Lax difference operator 



> 

•i-H 

X 

L = ]>> fc (x)A fc , (l.l 

where A is the translation operator A k f(x) = f(x + k)A k , or equivalently [A, x] = A. In 
these cases, the dynamics of the field Uk(x) are governed by a Lax equation 

or 

— = [V + (L n ),L], n>0, (1.2) 

where V + is a projection-like operator, whose detailed form is discussed in Section H~2l It is 
one of the hallmarks of the discrete hierarchy that the dynamical field Uk(x) only interacts 
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with itself in points that belong to the same (affine) integer lattice Z + x. Therefore, if one 
assumes that the Hamiltonians V + (L n ) do not have explicit x-dependence, one may ignore 
the fractional part of the space coordinate x, as it only labels isomorphic non-interacting 
substructures of the discrete hierarchy. In other words, it is legitimate to consider the space 
variable x to live on the set Z of integer numbers rather than the whole continuous space R. 
This trivial fact reduces the discrete hierarchy to an infinite, but countable, matrix problem. 
The matrix picture becomes even clearer if one writes the dynamical field Uk(x) as a matrix 
u i:j , i,j G Z: 

u k (x) = u x>x+k = u(x,x + k) or u {j = Uj-i(i) . (1.3) 

It is convenient to call Uk{x) = u X)X+ k a link from a lattice point x G Z to a lattice point 
x + k G Z, and to call the integer k G Z the (signed) length of the link. Translating the 
difference operator L into a matrix u = (My) turns out to be very fruitful, partially because 
this provides a better geometric and algebraic understanding. For illustrative purposes, 
consider the seminal bi-Poisson structure of Ref. OQ, here extended to both positive and 
negative link lengths 

{u n (x),u m (y)} K1 = 
{u n (x),u m (y)} K2 = 



and where e denotes the sign function. Notice that the Kupershmidt Poisson brackets {•, - }ki 
and {•, - }k2 are non-local in the x and y coordinates. It is quite elaborate to verify by brute 
force that {■, - }ki and {■, - }k2 are a bi-Poisson structure, not to mention identifying allowed 
deformations of this. We shall see that the non-locality is an artifact of the difference operator 
language, and that the matrix picture, combined with classical i?-matrix theory, provides an 
elegant and effective formalism to deal with such bi-Poisson structures. 

The outline of the paper is as follows. We continue this Section ^ with describing the con- 
nection between infinite matrices and difference operators. This provides a useful connection 
between the Lax formulations of the two pictures. In Section |21 we revisit classical i?-matrix 
theory for the matrix formulation, where in particular the conditions for integrability are 
given in detail. We provide local and non-local classes of solutions to these conditions, that 
are useful for model building. In Section El the Lax equations of motion are recast as the 
Hamiltonian equations of motion by finding viable Poisson brackets, that gives rise to bi- 
Poisson structures. It is natural in this connection to briefly review the three -R-theoretic 
Poisson brackets and the Suris quadratic bracket construction [2 . As a new result, we show 
that the Suris bracket {-, -}s2 decomposes into the second -R-theoretic bracket {-, and 
a piece {-, - }n2, that plays no role for on-shell dynamics. Our main results for the bracket 
structure are given in eq. 

flEZP - CEZfl and eq. (!3.3(-iaj) - (l3~3Kl) . In Appendix El we translate 
2 
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n, m G Z 



our Poisson bracket results into the notation of Kupershmidt for comparison. From a 
i?-theoretic perspective, it is remarkable that already the original discrete hierarchy of Ku- 
pershmidt induces the generalized construction of Suris. Finally in Appendix [BJ we give a 
2x2 dimensional example that illustrates aspects of the Suris theory. 

1.1 An Algebra Isomorphism 

We consider an algebra isomorphism from the associative algebra A of infinite-dimensional 
matrices (to be specified below) to a certain sub-algebra Aa of difference operators, 

A 3 u h-> L = u x , x+k A k G Aa ■ (1-5) 

fcez 

The algebra isomorphism maps matrix multiplication into composition of difference opera- 
tors: 

J2(uv) x , x+k & k = (J2 u x , x+kl A kl ) (J2 v x , x+k2 A k2 ) , u,v e A . (1.6) 

k&L feieZ k 2 & 

Geometrically, it is a useful fact that the total length k = k\ + k 2 of links is preserved under 
the matrix multiplication/composition. 

As an easy application, the matrix description provides a pictorial understanding for the 
coefficient functions {L n )k (x) of the higher powers of the Lax operator 

L " = E E u x , x+kl A kl ---u x , x+kn A k - = ]T(L") fe 0r)A fc , n>0. (1.7) 

fcez ki,...,k n ez kei 

They are 

n / 'i — 1 r \ 

(L n ) k (x)= Yl H u U + J^k^x + ^h = (u n ) XtX+k , (1.8) 

fci,...,fc„GZ r=l \ i=l i=l / 

k i -\- . . . -\~ k n = k 

or in words, the fc'th coefficient of the n'th power of the Lax operator corresponds to n 
consecutive links with the two free ends a distance k apart. Therefore, 

L n = ^T^WA*, n>0. (1.9) 

kez 

A difference operator, that has all the x's appearing to the left of all A's, is called a normal 
ordered difference operator. 3 

Let us now define the matrix algebra A itself. For simplicity, let A be the algebra of matrices 
with only finitely many non-zero entries 

A .= {it = (uij^Uij = for |i| + \j\ > 0} . (1.10) 
3 Ref. pQ uses anti-normal ordering, with all the A's appearing to the left of all a:'s. 
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This choice ensures that the matrix multiplication and the matrix trace (tr) are well-defined 
operations. 4 

1.2 Lax Formulation 

In this Section, we return to the Lax equation ()1.2j) and implement the corresponding matrix 
formulation, which appears naturally in the theory of classical P-matrices (cf. the next 
Section I2J). The Lax equation for L can be written in several ways 

BT 1 

— = [P + (L% L] = [L, V4L n )) = -[lZ(L n ), L] , L e A A , n>0, (1.11) 

for the operators V± and 1Z related through 

V+ + V- = 1 and K = V+ - V- . (1.12) 

Although we are going to consider different examples of the operator triples (V+, V-,7l), we 
will always assume that the three operators V+, V- and 1Z are interlocked via the above two 
relations (jl,12|) . Hence it is always enough to specify one of them. As an example, consider 
Hamiltonians of the form 

V + (L n ) := (L n )>o + ^(L n )o , n > , (1.13) 

where v is a constant parameter which is related to a choice of operator ordering pre- 
scription. The choice v — 1 leads to the standard hierarchy with the Hamiltonians given 
by (L n )>o, while v — — 1 leads to the so-called modified hierarchy with the Hamiltonians 

(L n )>i-~ 

To implement the matrix program, one seeks the matrix counterparts R, P± : A —>■ A of the 
operators TZ,V± : A a — > A a- They satisfy 

P + + P_ = 1 and R = P+ - P- (1.14) 

as well. Inspired by eq. (jl.9j) . we claim that the sought-for identification is provided by 

S (1.15) 
V±{L n ) = ^(P ± (^)),,, +fc A fc , n > . 

fcez 

It may look discouraging that (for instance) the P + operator does not act on the A-part at 
all, but only on the tt-part, as one usually counts the power k in the A fc -factor to determine 

4 If one forgets the associative matrix multiplication structure, but keeps the Lie commutator operation, 
A is often referred to as the Lie algebra gl(oo) in the mathematical literature. Also we stress that the above 
A has no identity matrix and no invertible matrices in this infinite dimensional case. At a few places in the 
paper we refer to invertible matrices and to make rigorous sense of this, A should be chosen as a Banach 
algebra with an algebra norm. It is out of scope to provide details here. 
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the action of V + (cf. eq. (|1.13j) ). However, one should recall that k is also available in the 
w-part as a link length. 

One can now lift the Lax equation (jl.llj) to the corresponding matrix algebra A: 

J£ = [P+(u n ),u] = KP-K)] = \[R{u n ),u], ueA, n>0. (1.16) 

In fact, the equivalence of the Lax eq. (jl.lljl for the difference operator L G A a arid the 
Lax eq. (jl. 16)1 for the matrix u G A follows straightforwardly from the algebra isomorphism 
()1.6|) and the prescription of eq. (|1.15|) . 

Let us work out the corresponding matrix maps R, P± : A — > A of the example eq. ()1.13|) 
given above. To this end, we need to introduce some notation. Let G A denote an 
elementary matrix such that (ey)^ = 8^8 jj. As is well-known, the e^-'s constitute a standard 
basis for the matrix algebra A, and hence a generic algebra element u can be decomposed as 
u = J2i jez u ij e ij e A. Similarly, one only have to determine the linear maps R, P± : A — » .A 
on the basis e^. The conversion prescription of eq. (|1.15jl is satisfied if one let 

R{eij) = e v (j-i) e i:j , 
P+faj) = 9 u {j-i) e l3 , (1.17) 
P-{eij) = 0- u (i-j) tij , 

where we have defined a sign function £j/i^x^ as 

{1 for x > 
1/ for x = (1.18) 
— 1 for x < , 

and a corresponding step function 

{1 for x > 
for x = (1.19) 
for x < . 

2 i?-Matrix Formalism 



The classical i?-matrix theory provides a universal method for constructing three compatible 
Poisson structures and infinitely many commuting charges for a wide class of integrable 
models jl]. 5 A classical i?-matrix is by definition a linear map R : A — » A such that the R 
bracket 

[u,v] R := i[i2(u),i;] + i[u,i2(i;)] (2.1a) 

=JP+(u),P + 00] - [P_(«),P.»] (2.1b) 

5 See also e.g. El El 13 El Ej • In HO) an extension of the i?-formalism to the fermionic Toda model is 
given. 
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is a Lie-bracket, i.e. it satisfies the Jacobi identity. A sufficient condition for the Jacobi 
identity is provided by the modified Yang-Baxter equation YB a (it, v) = 0, where the modified 
Yang-Baxter operator is given by 

YB a (u, v) := [R(u),R(v)]-2R[u,v] R + a[u,v] . (2.2) 

As we already have seen in the previous Section, it is convenient to define projection-like 
operators P± : A — > A: 

P ± := i(l±i?) , R = P+ ~ P- ■ (2.3) 

We emphasize that in general Im(P + ) and Im(P_) do not form a direct sum, and P± are not 
necessarily idempotent operators. 

It is instructive to see how the integrable model arises. Consider an abelian subalgebra 
Ao ^ A. (Usually we simply consider the infinite hierarchy A® D { u n \ n— 1, 2, . . .} generated 
by a single algebra element u E A.) The dynamical flow 

5 v u = [P+(v),u] = [u,P-(v)] = -[R(v),u] , u,veA , (2.4) 

is generated by a Hamiltonian P + (v). After some straightforward algebra, the commutator 
of two flows reads 

[6 w ,S v ]u = ^[N(v,w),u] , u,v,weAo, (2.5) 

where we have defined the Nijenhuis tensor [TI] 

-N(u,v) = N + {u,v) + N^{u,v) (2.6a) 

= [P±(u),P±(v)]^P±[u,v] R , (2.6b) 

and the chiral Nijenhuis tensors 

N ± (u,v) := P T [P ± (u),P ± (v)} . (2.7) 

The Nijenhuis tensor is equal to the modified Yang-Baxter operator 

N{u,v) = [u,v} + [R(u),R{v)}-2R[u,v} R = YB^v) . (2.8) 

The flows S n , n > 0, commute for an integrable system. From eq. (J2.5|) . a sufficient integra- 
bility condition is provided by the modified Yang-Baxter equation 

YBi(«,v) = N(u,v) = , (2.9) 

which our examples in Sections 12. lH2.2l will satisfy. More generally, integrability is guaranteed 
if there exists a linear operator B : A —>■ A, such that 

N(u,v) = B[u,v] , (2.10) 
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as can easily be checked from eq. (|2.5j) . When we assume a vanishing Nijenhuis tensor 
N = 0, it follows from eq. (|2.6b[) . that the ±P± operators are Lie-algebra homomorphisms 
(A, [•, •]#) — > (.4., [•, •]), and in particular that the images Im(P ± ) are two Lie sub-algebras: 

[Im(P ± ),Im(P ± )] C Im(P ± ) . (2.11) 

The vanishing of the chiral Nijenhuis tensors 

N + = and AL = (2.12) 

implies the vanishing of the Nijenhuis tensor N = YBi = (cf. eq. (|2.6aj) ). The opposite 
statement is not true. The local and non-local examples in the next Sections I2.1H2.2I will 
meet the stronger condition (|2.12j) . 

Both the vanishing Nijenhuis tensor condition ()2.9|) and the vanishing chiral Nijenhuis tensor 
condition ()2.12|) are stable under conjugation of the P-matrix R'(u) = aR(a~ 1 ua)a~ 1 with 
an invertible algebra element a. The conjugation procedure can be used to generate new 
P-solutions from old P-solutions, although we will not pursuit this here (cf. footnote 0}. 

2.1 A Class of Local R-matrix Solutions 

Here we propose a class of local solutions P, P± : A — » A to the condition eq. ()2.12j) that is 
parametrized by an arbitrary function v : Z — > C and that generalizes (|1.17j) . It is given by 

P(ey) = £ v (j,i) e i:j 

P+M = Q„(j,i) eij (2.13) 
P-(e lj ) = Q- U (i,j) , 

where we have defined a generalized sign function (involving now v = u(x) being a local 
function) 

{1 for x > y 

v{x) for x = y (2-14) 
— 1 for x < y 

and a corresponding step function 

1 for x > y 

e„(x,y) := l[l + £ v (x,y)} = ii±|M f OIX = y (2.15) 



2' 



2 

for x < y . 



The choice of v : Z — > C is related to a (x-local) choice of operator ordering prescription. Here 
we refer to a P-solution as being local if the standard basis diagonalizes the P-matrix. 
(We emphasize that a local solution usually becomes non-local in terms of the difference 
operator fields Uk(x).) We prove in the next Section f2.2l that the chiral Nijenhuis tensors 
vanish (cf. eq. (|2.12j) ). so that P satisfies the modified Yang-Baxter equation YBx(P) = 0. 



7 



The operators P + and P_ from eq. (|2.13|) "project" onto (weakly) upper or (weakly) lower 
triangular matrices, respectively, but they may share diagonal matrices unless u(x) = ±1 : 

Im(P+) n Im(P_) C {0} <S> Im(z/) C {±1}. (2.16) 

Thus in general, the sub-algebras Im(P + ) and Im(P_) do not form a direct sum. Similarly, 
the operators P + and P_ are idempotent (P± = P±) if and only if the P-matrix is an 
involution (P 2 = 1), which holds precisely when v{x) = ±1, as expressed by 

4 P+ P_ = 1 - P 2 = [l-v 2 ) 5 , (2.17) 

and where 5 : A — > A projects onto the diagonal matrices 

5( eij ) ■= Sij etj , 5 2 = 5 , P 5 = v 5 = 5 R . (2.18) 

2.2 A Class of Non-Local i?-Matrix Solutions 

There is a non-local generalization [3] of the solutions in ()2.13|) that reads as 

Pf(Ci i j) ^(j ^) 6jj + 2^''' ^ ^ ^i,m Cmm 

P - (^-ij) $(^ j) Cjj 2 < ^*'"' ^ ^ ^mm 

mgZ 

The diagonal case z/jj = i^dJy corresponds to the previous local solution. We claim that the 
non-local P-matrix possesses vanishing chiral Nijenhuis tensors (cf. eq. (|2.12Jl ). To prove 
this, it is enough to consider iV±(ey,ejy) for two basis elements and e^. P + and P 
"project" onto links with weakly positive and weakly negative link length, respectively. The 
Lie-bracket preserves the total link length. So to give a non-zero contribution to N±(eij, e^i) 
both entries and e^i have to be zero-length links. But the zero-length links are nothing 
but the diagonal matrices and those commute trivially. 

In the non-local case, the sub-algebras Im(P + ) and Im(P_) form a direct sum if and only if 
the matrix Uij is an involution: 

Im(P+) n Im(P_) C {0} v 2 = 1 . (2.20) 

2.3 The i?-Bracket 

For the local and non-local solutions (J2.13|) and (J2.19j) . the P-matrix P = + R^ is a 
linear function of u, where the superscript "(0)" and "(1)" refer to the power of v. The 
P-bracket [• , - }r inherits this linear //-dependence, and can be split accordingly 

[•>•]* = h'ff + h-ff (2-21) 

8 



(2.19) 



into two mutually compatible Lie-brackets [ • , ■ ]^ and [ ■ , • ]^. By definition, the i?-bracket 
is a Lie pencil in v. We have 

[eij, e M ) R = - [£ v (j, i) + £ v (l, k)} (S jik e a - 5 itl e kj ) ^ 
= [9„(j,i) 9v{l,k) - ®- u {k,l)) (8 j>k e u - 8 it i e kj ) 

for the local .R-matrix (j2.13j) . This is a gl(oo) Lie algebra [e^, e k i] = 5j ik eu — &i,iZkj with a 
2-cocycle-like prefactor. For the non-local R- matrix (j2.19|) . only the first-order contribution 
in v is changed. It reads 

[ e ij, e ki\ ( R = 7^$i,j {vj,k ~ v it i) e M + -8 kj i (u ktj - vij) . (2.23) 

It is a curious fact that links e mm of zero- length can never be produced in a i?-bracket [ • , • }r 
of the local or non-local type (cf. eqs. ()2.13|) and (j2.19|0 . Specifically, 

ti(e mm [u } v] R ) = . (2.24) 

In contrast, the same does not hold for the standard gl(oo) Lie-bracket [•, •], where for 
instance [e^, e 3 j] = eu — ejj yields two zero- length links if i ^ j. 

2.4 Equations of Motion and Time Evolution 

It is interesting to write out the equations of motion in coordinates. If we insert the local 
solutions (j2.13|) for R, P + and P_ into the Lax eq. (jl.l6j) . we get 



dt n ' 2 

n kez 



2 ^2 i) {u n ) ik u kj - u ik S u {j, k) {u n ) kj } (2.25a) 

fegZ 

J2 [©.,(*, i) (u n ) ik u kj - u lk Q u (j, k) (u n ) kj ] (2.25b) 

k&Z 

[ u « ®-v( k i3) Wkj ~ ®-v{h k) (u n ) ik u kj ] , (2.25c) 



kez 

respectively. The non-local generalization of eq. ()2.25aj) reads 
dun 



dt n 2 ^ 



]- Y [e(k-i) {u n ) ik u k j - u ik z{j-k) {u n ) kj ] 



+ ^Uij ^2(vk,i - v kyj )(u n ) kk . 
kez 

The non-local generalizations of eqs. (j2.25b|) and ()2.25c|) are similar. 



(2.26) 
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2.5 Lattice Truncations 



We next address the question whether we can constrain the dynamical fields Uk(x) without 
violating the equations of motion? As is easily seen in the matrix formalism, it is consistent 
with the equations of motion ()2.26|) to deploy "rectangular" type of truncations of the matrix 
algebra A to a sub-algebra 

Ai := {u = (u i:j ) e A | Uij / i,j e 1} (2.27) 

for some index-set /CZ. From a difference operator perspective, it is natural to consider 
"diagonal" type of truncations. Here, we discuss two "diagonal" type of truncations that are 
often used in applications: 

• Truncation of the link length from below: It is consistent with the equations of motion 
to consider a truncated model with 

\fk < N : u k {x) = . (2.28) 

To prove this, we note that the first (second) term on the rhs. of (12.25b)) has an index 
variable k > i (k < j), so that the left hand side duij/dt n depends at least linearly on 
a link u k j (uik) with a signed length j — k<j—i(k — i<j — i). So if the field (and 
its fellow fields with less or equal link length) are annihilated at some time (t\, ti, ■ ■ ■), 
the equations of motions (|2.25bJ) cannot undo that for other times. 

• Truncation of the link length from above: A similar examination of eq. (|2.25c|) shows 
that there is also a consistent truncation from above: 

\/k> M : u k (x) = 0. (2.29) 

The two truncation schemes are also consistent with the non-local solution ()2.26|) . because 
when considering the left hand side dv,ij/dt n , the non-local terms are always hidden behind 
at least one power of u^. By invoking both of the above truncation schemes, we get models 
with only a finite number of different fields . . . ,un(x) with link lengths between M 

and N; all located inside an infinite universal enveloping construction. This fact renders the 
discrete hierarchy highly accessible for applications. 



3 Poisson Brackets 



Before we proceed with constructing Poisson brackets, we need to introduce a few standard 
notions to fix the notation. A non-degenerate bilinear form 

(u,v) = tr(-ut>) = (v,u) (3.1) 

is inherited from the matrix trace (tr). Note that a fo-linear form (•,•), in contrast to a 
sesgm-linear form, has no internal transposition (or Hermitian conjugate for that matter). 
This is mainly to ensure that the bilinear form is invariant/associative: 

(u, [v,w]) = ([u,v] , w) . (3.2) 
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The non-degenerate bilinear form gives rise to an identification of the algebra A with the 
set A* of linear functionals on A. For a linear operator R : A — > A, the dual operator 
R* : A* — > -4* becomes identified with the transposed operator 

(v,R(u)) = (R*(v),u) . (3.3) 

One may always decompose an operator R in symmetric and skew-symmetric parts 

_ R i R* . . 

R± ■= —j— • (3-4) 

Let us note for later that the non-local i?-solution (J2.19)) from Section 12.21 decomposes as 

R -( e ij) = e {j-i) e ij + 7,Ki ^ V[iM e mm (3.5) 

R+ifiij) = 2^ i '~' ^ ^ ^{i,m}&mm ■ (3-6) 

Notice that the skewsymmetric part i?_ in this case is again a i?-matrix, with vanishing 
chiral Nijenhuis tensors (cf. eq. ()2.12|) ). 

The adjoint action ad(w) : A — > A is defined as ad(u)v := [u,v]. Because of the invari- 
ant/associative property of the bilinear form, the coadjoint action (from right) ad*(w) : 
A* — > A* is identified with minus the adjoint action 

(&d*(u)v,w) = {v,ad(u)w) = (v,[u,w]) = — ([u, v],w) = — (ad(u)v, w) . (3.7) 

The gradient V/ of a function f = f(u) on the dual space u G A* can be defined implicitly 
via the infinitesimal variational formula 

6f = (6u,Vf) . (3.8) 

Explicitly, the gradient is 

dun 



7« 



*J6 

Notice the transposition of indices in the above formula. 
3.1 Conserved Charges 

As is well-known, a hallmark of an integrable system is an infinity of conserved charges. In 
the discrete hierarchy, the charge densities are defined as 

h n (x) = -(L n )o(x) = -(u n ) x , x , n>0, (3.10) 

n n 
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and the charges themselves are defined as 

H n = J2h n (x) = -tr(w n ) , n>0, (3.11) 
xez n 

where the trace in the last equation can be thought of as the sum of all closed loops build 

out of n consecutive links. It follows directly from the Lax eq. 16|) and from the cyclicity 
of the matrix trace tr[tt, v] = 0, that the charges H n , n > 0, are conserved in time. Also, one 
observes easily that 

(L n ) k (x) = (u n ) x , x+k = , (3.12) 

or, equivalently, 

VH n+l = u n . (3.13) 

The above facts are slightly more elaborate to establish purely within the difference operator 
approach (cf. for instance Theorem III. 2. 7 of Ref. pQ). 

One may rewrite |H |Hj the Lax eq. ()1.16|) as a standard Hamiltonian equation on the dual 
space A* 

^ = -[RVH n+1 ,u] = --ad*(RVH n+1 )u, u e A* . (3.14) 

It is useful to seek for Poisson bracket structures {■, -} p , p = 1, 2, . . . , n + 1, such that the 
time evolutions can be reproduced as Hamiltonian flows with the conserved charges H n+2 - p 
acting as generators 

{f(u),H n+2 . p } p = «P (^V/) °P \([R{rr)M, V/} . (3.15) 

In other words, find Poisson brackets such that the Lax eq. (jl,16|) can be written as the 
Hamiltonian eq. (J3.15)) . This will ensure the Lenard relations 

{•,H n+1 } p = {-,H n } p+1 , n,p>0. (3.16) 

Note that the H n 's play a double role as both Hamiltonian generators and as conserved 
quantities. This forces them to mutually "commute" in a Poisson sense: 

{H n ,H m } p = 0, n,m,p>0. (3.17) 

For the local and the non-local solution discussed in Sections 12.1112.21 each Poisson bracket 
{-, - } p will be a linear function of ^-parameter 

{;-} P = {;-} P ° ] + {•,•}« (3-18) 

of two mutually compatible Poisson brackets {•, - } p °^ and {•, -}p ■ This is sometimes referred 
to Poisson pencil in v. 
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3.2 1st Poisson Bracket 



We now derive the first Poisson bracket from the i?-matrix formalism. From eq. (|H.15jl and 
the Lax eq. (jl.lfij) . we get 

= l([R(un,u],Vf)+l([u n ,u],R(Vf)) (3ig) 

= l{u,[Vf,R(u n )]) + ~{u,[R(Vf) 7 u n }) 
= (u,[Vf,u n ] R ) = (u,[Vf,VH n+1 ) R ) , 

where the second term (which is identically zero) was added to achieve the required skewsym- 
metry. We immediately recognize the first .R-theoretic Poisson structure 0] 

{f,g}m ■= (u,\yf,Vg] R ) ■ (3.20) 
For the local i?-matrix solution (|2.13j) . this leads to 

{uij,Uki}m = (u, [eji,ei k } R ) = ^ [S u (i,j) +£ v (k,l)] (6 itl u kj - 5 jtk u u ) . (3.21) 

The local bracket {•, -} R i has a pictorial interpretation as a concatenation of two "incoming" 
links with a 2-cocycle-like prefactor. For a non-local v, the first-order contribution in v is 

{uij.Uki}^ = -Sij {v it i - is j)k ) u ki + ~8 kj i (ui ti - v k j) Uij . (3.22) 

The bracket operation {•, -}^] preserves at least one of the "incoming" links «y and u k i (up 
to an overall factor). To have a non- vanishing "outgoing" link, say u^, the other "incoming" 
link u k i should have zero link length, and it should "interact at a distance" with an endpoint 
of the first "incoming" link. The interaction at a distance is mediated through a non- 
vanishing, non-local matrix element v n ^ m . In the difference operator language, the local 
bracket transforms into 

{u n (x),u m (y)} m = - [s(n) - v x 5 nfi + e(m) - v y 5 mfi ] ^ 

^ [^n+m^) &x+n,y "n+m(2/) ^x,j/+r?i] • 

3.3 Higher Poisson Brackets 

We now generalize the method of Section l3~21 to the higher brackets. Postponing the question 
of the Jacobi identity, let us tentatively write down skewsymmetric candidates for the (p + 
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l)th bracket structure {•, -} p+ i tr , labeled by an integer r = 0,1, ... ,p. From the eq. ()3.15|) 
and the Lax eq. we get 

{f, H n+1 _ p } p+ i r = — j- = (— — , V/) 
ot n dt n 

= h[R(u n )M,Vf) + l([u n - p ,u],R(u r (V/) u*- r )} 

2 2 (324) 

= -(u, [Vf,R(u r u n ~ p tiT*)]> + 2<«. M«' (W) vF- r ),u n -v]) 

= ~ (it, [V/,i2(u r (Vtf n+1 _ p ) rO]> + [J2(t* p (V/) O.VfWi-p]) , 

where the second term (which is identically zero) was added to achieve the required skewsym- 
metry. In this way we obtain the r'th candidate for the (p+ l)'th i?-matrix Poisson structure 



{/,JU := 2<«' [W,i2(« P (Vp) u p - r )]) + ~(u, [R(u r (V/) tT), Vs]> . (3.25) 

The eq. f!3.24j) is an inhomogeneous linear equation in the Poisson structure, with source 
terms generated from the Lax eq. (jl.lfij) . We may not have properly identified possible 

( H) 

homogeneous Poisson bracket parts {-, -} p+ i that commute with all the charges H n , n > 0. 
Besides these homogeneous contributions, the potential bracket candidates are convex linear 
combinations of the basis brackets {■, -} p+ i tr , r = 0,1, ... ,p. Again we stress that most of 
brackets are going to be discarded, as they will not satisfy the Jacobi identity. 



3.4 2nd Poisson Bracket 



As we saw in the last Section, the potential bracket candidates for a quadratic bracket 
include the convex linear combinations of the two basis elements {•, -}2,o and {-, -}2,i- The 
2nd -R-theoretic Poisson structure |3] turns out to be the symmetric average 

U,9}r2 ■= \{f,9h,o + \{f,9h,x = \(u,[Vf,R{u,Vg} + ])-(f^g), (3.26) 

where {u, v} + := uv + vu. If the Jacobi identity holds, one may show that {•, -}# 2 is always 
compatible with the first i?-theoretic bracket {-, -}_ri. This follows by shifting u — > u + Al 
in eq. (|3.2f)jl . because the shifted 2nd Poisson bracket 

{f,g}m{u+Xl) = {f,g}m(u) + X{f,g} m (u) (3.27) 

can be re-interpreted as a Poisson pencil between the two brackets. (It is enough to let / and 
g be linear functions of u, so that V/ and Vg are w-independent, which simplifies the above 
argument.) Moreover, one may prove |3] that sufficient conditions for the Jacobi identity 
for the {-,-}r2 bracket are, that R and R_ satisfy the modified Yang-Baxter equations 
YB a (R) = and YB a (i?_) = with the same parameter a. This indeed is the case for the 
local and non-local solutions (cf. eqs. ()2.13|) and (|2.19|) ). 



14 



A generalization of the quadratic bracket is due to Suris 0. He defines a bracket 

2{f,g} S2 := (A 1 ((Vf)u),(Vg)u)-(A 2 (uVf),uVg) 

+ (Si(uV/),(V0) u) - (S 2 ((Vf) u),uVg) , 1 " J 

where Ai, A 2 , Si, and 5*2 are linear maps A — > A satisfying A* = —Ai and S 1 * = S 2 . We 
assume everywhere in Section^ that 

R = A 1 + S 1 = A 2 + S 2 , (3.29) 

and that both A\ and A 2 satisfy the modified Yang-Baxter equation YB a (Aj) = 0, i = 1, 2. 

With the above assumptions one can show that the two Suris conditions 

2S 1 [u,v]a 2 = [S^tSiiv)] ^ 3o ^ 
2S 2 [u,v] Al = [S 2 (u),S 2 (v)} 

are sufficient for the Jacobi identity to hold. Also, they imply the modified Yang-Baxter 
equation YB Q (_R) = 0, and that {•, -}ri and {•, -}s2 are compatible Poisson brackets. 

Note that the opposite does not hold, i.e. that YB a (_R) = does not necessarily imply 
the Suris conditions ()3.30|) . We give a counterexample in Appendix |Bj Also, the 2nd 
i?-theoretic Poisson structure {•, - }r 2 can be seen as a special case of the Suris construc- 
tion if one let Ai = A 2 = R_ and Si = S 2 = R + , because in this case the Suris con- 
dition 2R + [u,v]n_ = [R + (u), R + (v)] does follow from the modified Yang-Baxter equations 
YB a (R) = and YB a (i?_) = 0. 

It is known that a compatible quadratic Poisson structure for the discrete hierarchy is not 
unique [2], although a full classification of ambiguities is still an open problem. Here, we 
give a family of solutions that can be described using the quadratic Suris bracket. To this 
end, define a skewsymmetric linear map Q : A — » A 

n := = = -Q* , (3.31) 

where we used eq. (13.29(1 in the second equality. One may decompose the Suris variables 
entirely in terms of R and Q: 

A 1 = R_ + tt , A 2 = - tt , S 1 = R + ~n , S 2 = R + + Vl , (3.32) 

as well as the Suris bracket itself 

{v>S2 = {;-}m + {;-}a2 , (3.33) 

where 

{f,g}n2 ■= ^0[ Ul V/],[u,V 5 ]) . (3.34) 

The structure {•, -}q 2 is not required to satisfy the Jacobi identity, and hence it is not 
necessarily a Poisson bracket, although this turns out to be the case for our example below. 
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A sufficient condition for this to happen is given by the Yang-Baxter equation YB (fi) = 0. 
In general, the structure {•, -}fi2 does not contribute to the Hamiltonian eq. ()3.15j) . because 
{H n , ■ }q 2 = 0, n > 0. Hence, even for the more general Suris bracket {•, -}s2, the dynamics 
is governed by the {-, - }r2 bracket alone. 

We claim that the non-local i?-solution ([2.19)1 . together with the choice 

(3.35) 

for some skewsymmetric matrix Uij = meets all the conditions of the Suris construc- 

tion. The proofs are very similar to the discussion given in Section 12.21 First of all, both 
A\ and A2 are of the same form as the non-local i?-solution ()2.19|) . and therefore they too 
have vanishing chiral Nijenhuis tensors N±(Ai) = 0, and hence YBi(Aj) = 0. Secondly, both 
sides of the Suris conditions (|3.3(Jj) vanish. For instance, the lhs. is of the form Si(w), where 
w = 2[u,v]a - Because of the special form of the two Si maps, i — 1,2, only diagonal parts 
of w could potentially contribute. On the other hand, diagonal parts of [u, v]^ do not exist 
according to eq. (j2.24|) . So the lhs. is zero. The rhs. is zero, because both Si(u) and Si(v) 
are diagonal matrices, and hence commute. 

Let us write down the Suris quadratic bracket {•, -}s2 = - fsl + {'> '}s2 * n detail 

{uij,Uki}sl = \ [e(k-i) +e(l-j)]u u u kj , (3.36a) 

{uij, u k iYsl = Ld ijM Uij u k i , (3.36b) 

where 

Vij,kl '■= ~Ay{i,l} ~ u {j,k} + V[k,i] + + + U it k+ U k j) = — LOkl,ij ■ (3.37) 

For simplicity, we have collected all the t^j-terms inside the {•, - j^-part. (Strictly speaking, 
this represents a minor abuse of notation, because Ui j does not need to be first order in 
v.) Figuratively speaking, the Suris bracket {•, -}s 2 consists of two parts {■, and {■, -j-^ 
that represent elastic and inelastic scattering of two "incoming" links, respectively. In other 
words, the first order bracket {•, preserves the two "incoming" links Uij and Uki, while 
the two "incoming" links exchange a pair of endpoints in the zero order bracket {•, -}gl- 

If we restrict the ^-contribution to be the form Uij = uu^j = —cjj-i, the quadratic bracket 
reads 

{u n (x),u m (y)} { gl = i £ ( k ) + £{k+m-n)]u n - k (y) u m+k (x) 5 x+k , y 

fcez 

{u n (x),u m (y)}gl = \ u n{x) u m (y) y^jy x h,-m - Vy h,n (3.38) 

~ ^k+m-n + ^k+m ~ ^k + ^k-n] &x+k,y 

in the difference operator language. 
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3.5 3rd Poisson Bracket 



The potential bracket candidates for a cubic Poisson bracket are convex linear combinations 
of the three basis elements {•, -^o, {■, -}3,i and {-, -}3,2- The 3rd -R-theoretic bracket IE] 
turns out to be given entirely by the candidate {-, -^i.: 

{f,g}m := {f,9h,i = \(u, [Vf,R{u (Vg) u))) + ^{u, [R(u (V/) u), Vg}) , (3.39) 

One may show E] that the three -R-theoretic brackets {-,-}m, {', 4 }r2 an d { 4 , 4 }m are 
compatible Poisson structures if both R and _R_ satisfy the modified Yang-Baxter equation 
YB a (R) = and YB Q (RL) = with the same parameter a. In general, the third bracket 
{-, - }r 3 is not compatible with the Suris quadratic bracket {-,-}s2- 

We derive 

{uij,u k i} R3 = - ^[£„(m, i) + £ V {1, m)] 

™ (3.40) 
- 7, ^2[ £ »tii m ) + £ v{m, k)) Uim U m l Uj~j 

m 

for the non-local .R-solution ()2.19|) . Its local first order v terms are 

{uij,u M }xl = ^Uij [u ki Vi uu -u kj Uj uji] - <-> (k,l)} , (3.41) 

while the non-local first order v terms read 

{Uij,U k l}xl = ^Ujj ^ Ukm (v m ,i - Vm,j) Uml - [(h j) <-> (k, I)] ■ (3.42) 

There is an interesting duality between the 1st and the 3rd bracket, which (formally) facil- 
itates the proof of the Jacobi identity for the third bracket. Following Oevel and Ragnisco 
one notices that matrix inversion u \— > u^ 1 maps the first and third bracket into each 
other up to an overall minus sign. In detail, consider linear functionals f(u) = (a,u) and 
g(u) = (b,u) for some constant algebra elements a,b e A. Then for invertible w's 

Vf = a, Vg = b, (3.43) 

V(/(m" 1 )) = -u^au' 1 and V(^(w -1 )) = - u^bu' 1 . (3.44) 

It follows from the definitions (|3.20j) and (|3.39|) that 

{f(u- 1 ),g(u- 1 )} m (u)= ^w-W 1 ], /?(&)> - (a<-6) 

= ±([a,u- l ],R(b)) - (o«-6) ^ 
= (t*- 1 ) ■ 
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This provides a proof of the Jacobi identity within the matrix group of invertible matrices 
(cf. footnote EJ). Similarly, one notices that the 2nd bracket is self-dual under u i— > u~ x up 
to an overall minus sign 

{f(u~ 1 ),g(u~ 1 )}i a (u)= ^{[u,u~ l au~ l ] , i2{u, u^&iT" 1 }*) - (a 6) 

= i^tr 1 ],^,*- 1 }^ - (a<-&) ^ 3 - 46 ) 
= -{f,g} R 2(u- 1 ) . 

Interestingly, the structure {-, -}n2 is also self-dual under u i— > m^ 1 , but with an overall plus 
sign 

{/(n- 1 ),^( M - 1 )} n2 H= ^(fi^tt-W],^^- 1 ^- 1 ]) - (a <-&) 

= +{f,g}n2{u- 1 ) , 

so different parts of the Suris bracket {-,-}s2 has different transformation properties under 
duality. Of course, one may claim the Suris bracket is self-dual under u t— > u~ l if one 
simultaneously changes the sign of Q, or equivalently, one simultaneously exchanges A\ <-> A 2 
and Si <-> S 2 - 

Acknowledgment: This work has been partially supported by DOE grant DOE-ER-40173. 



A The Kupershmidt Bi-Poisson Structure 

Here we translate our results into the q k (x) fields used by Kupershmidt [T]. He uses an 
anti-normal ordered Lax operator of the form 

L = ^A- fc g fc (x) . (A.l) 

fcez 

Comparing with eq. (jl.ljl . one derives the translation formula 

q k (x) = u k+XtX = u- k (k+x) . (A. 2) 

We may facilitates the u k (x) <-> q k (x) translation of the Poisson structures by the following 
observation. First of all, from the matrix definition u k (x) = u XjX+k , one notices that a change 
of variables u k (x) <-> q k (x) corresponds to a transposition of the link matrix Uij <-> Uji. 
Secondly, notice that the {-,-}ri bracket (j3.21|) and (J3.22JI and the {•, -}s2 bracket (j3.36a|) 
and §MB 

are invariant under a transposition u — > u of the w-matrix combined with a 
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change of the sign v — > — v of the ^-matrix. Hence the u k (x) <-> q k (x) translation simply 
amounts to a change of the sign of v. The {•, -}ri bracket (|3.23j) and the {•, -}s2 bracket 
(I3.38j) become 

{q n (x), q m (y)}m = - [e(n) + u x 5 nfi + e(m) + v y S mfi ] 

^ \qn+mix) $x+n,y Qn+mi^V) $x,y+m\ i 

{q n (x), q m (y)}sl = 2 2 + £ ( k + m - n )} Qn-k{y) Qm+k{x) S x+ktV , ,^ ^ 

{q n (x), qm(y)}sl = ifinix) q m {y) J^t - ^ + v v 

kez 

— ^k+m-n + ^k+m ~ ^k + C^fc-n] $x+k,y , 77,, 771 £ Z . 

In the original model of Ref. pQ, the fields corresponding to positive link lengths are con- 
strained 

q-i(x) ~ 1 and V7c < -2 : q k (x) ~ , (A.4) 
so that the Lax operator reads 

L = A + ^A- fc g fc (x) . (A.5) 

fc>0 

The constraints have to be consistent with the equations of motion, written either as Hamil- 
tonian equations (j3.15j) or as Lax equations f)l. 16|) - with or without use of Poisson brackets, 
respectively. Previously in Section 1231 we saw that the constraints Vk < —2 : q k (x) ~ are 
consistent with the Lax formulation. Also, it is easy to check from the bi-Poisson structure 
eq. (jA.3|) that the constraints \/k < —2 : q k (x) ~ decouple from the theory in the Hamil- 
tonian sense, i.e. that the Hamiltonian vectorfields {q k {x), • } — vanish for both brackets 
when k < —2. On the other hand, the constraint q~i(x) ~ 1 induces non-trivial conditions 
on the model. From the Lax eq. (Jl.llj) using V- and the expansion eq. (J1.7)) . one derives 

= ^Zlq-xix) [(L») (x) - (L») (x-1)) , (A.6) 

so consistency requires v — 1. Moreover, in the Hamiltonian formulation, where one imposes 
that the field g_i(x) = l "Poisson commutes" with the other fields q n {x), n > 0, one is lead 
to the choice v — 1 and u k = kc — e(k) = —uj- k with some immaterial constant c £ C. 
(Again we stress that the on-shell dynamics are not affected by the {■, -}q2 contributions.) 
With this choice, the brackets (|A.3|) simplify to 



{q n {x),q m {y)}Ki qn+m\X) &x+n,y <?n+m(Z/) &x,y+m j 
{q n {x),q m {y)} { K2 = + e(k+m-n)]q n _ k (y) q m+k (x) 5 x+k)V 



2 

kez 



{q n {x),q m {y)} { K2 = ~<ln(x) q m {y) Y^i £ ( k ") +z{k+m-n 



(A.7) 

j„(.t?) a^{ii\ > \fA t ) 4- fA k,-\-m, — n ) 
fcez 

- e(Hm+-) - e(&— ra— -)] 5 x+k>y , n,m > 
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which agree with the formula (III. 3. 4) and the formula (III.4.15a-c) in Ref. pQ. 



B Example: Mat 2x 2(C) 

Here we give a counterexample, that shows that YB a (R) = YB a (Ai) = Y~B a (A2) = 0, taken 
together with the relation R = Ai+S{, does not necessarily imply the Suris conditions ()3.3()j) . 

Consider the 4-dimensional associative algebra A = Mat2 X 2(C) = C 4 . A convenient basis 
is given by the 2x2 unit-matrix cr 4 = 1 and the three Pauli <jj matrices, i = 1,2,3, which 
satisfy the relation 

3 

GiGj = Sijl + i^eijkCTk , i,j = 1,2,3, (B.l) 

fc=i 

where is the 3-dimensional Levi-Civita symbol. A non-degenerate, associative/invariant 
bilinear form is inherited from the matrix trace (tr): 

ov) = tx{a lx a v ) = 2S^ U , //, v = 1, 2, 3, 4 . (B.2) 

We now search for i?-matrix solutions to the modified Yang-Baxter equation YB a (R) = 0, 
where a G C is a given fixed complex number. Let us consider a linear injective map $ 

C 3 3 f=(r u r 2 ,r 3 ) £ Re End(A) , (B.3) 
that maps a complex rotation vector f into its rotation matrix R 

3 

R(ai) := % ^ ejjkrjQ-k , i = 1,2,3, 

R(l) := . 

In other words, R rotates the basis of (7j matrices, i = 1,2,3, around the rotation axis f. 
Note that R = —R* is skewsymmetric. The i?-bracket reads 

[<T U <Tj] R = a {i r j} , i,j = 1,2,3 , 

[i,-U = o. (R5) 

We claim that YB^.^(R) = 0, i.e. that 

2R[a IJi ,a u ] R = [R(a^),R(a v )]+ f ■ r [a^a v ] , /i, v = 1, 2, 3, 4 , (B.6) 

where f ■ f := Yli=i r i * s a "bilinear" length-square, i.e. without a complex conjugation. 
Technically, since the ovsector is trivial, the eq. (jB.6|) reduces to a "dual" Yang-Baxter 
identity 

3 

e ijk ([R(crj), Rjcrk)) - 2it! [<jj, o- k ] R + r • f [o-j, crfc]) =0, z = 1,2,3, (B.7) 
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which is easy to verify by direct calculation. 

Now let us apply this fact to a specific example. Define five rotation vectors 

r = (l,i,Va) , si = (1,^,0) = - s 2 . (B.8) 

and 

ai = r — si = (0,0, y/a) , a 2 = r — s 2 = (2,2i,y/a) . (B.9) 
The corresponding five rotation matrices are skewsymmetric 



Ai = $(a*) = - A* , 
Si = $(si) = $(-s 2 ) = - $(s 2 ) = - s 2 = s; , (B.10) 
i? = $(f) = $(ai + Si) = $(a i ) + $(s l ) = Ai + Si, i = 1,2, 

and they each satisfy a (modified) Yang-Baxter equation 

YB a (i?) = 0, YB a (Ai) = 0, YB (^) = 0, i = 1,2. (B.ll) 

Using YBo(Si) = 0, one may reduce the Suris operator 

i- 2S l [u,v] A2 -[S l (u),S l (v)} = 2S l [u,v] R (B.12) 

to only one term (cf. eq. (|3.3Uj) ). Therefore, the "dual" Suris condition simplifies to 

^3 3 

2 e ijkSi[crj,crk}R = i/^jS^y erg = - i f ■ 3 si )f 7^ 0, i = 1,2. (B.13) 
j,fc=i j=i 

Thus the Suris conditions ()3.30|) are not met, despite eq. (|B.11|) . 

It is however generally valid that the two Suris conditions eq. f)3.30|) . the modified Yang- 
Baxter equation YB a (R) = 0, taken together with the relation R = Ai + Si, imply the two 
modified Yang-Baxter equations YB Q (v4j) = for the skewsymmetric maps Ai, i — 1,2 (cf. 
Theorem 2 in Ref. [3]). 
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